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Abstract 

We introduce a comprehensive scheme to physically quantify both viscous and elastic rheological 
nonlinearities simultaneously, using an imposed large amplitude oscillatory shear (LAOS) strain. 
The new framework naturally lends a physical interpretation to commonly reported Fourier coef- 
ficients of the nonlinear stress response. Additionally, we address the ambiguities inherent in the 
standard definitions of viscoelastic moduli when extended into the nonlinear regime, and define 
new measures which reveal behavior that is obscured by conventional techniques. 
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Biopolymer networks [H [21 [3], wormlike micelles [1], colloidal gels [5], and metastable 
soft solids in general [B], can be classified as nonlinear viscoelastic materials and as such 
have been of interest to experimentalists for many decades (e.g. [7]). The biological and 
industrial processes associated with these materials often involve large deformations, yet 
standard methods of characterizing their nonlinear rheological properties rely on techniques 
designed for small strains. In this Letter, we develop a new and systematic framework for 
quantifying the nonlinear viscoelastic response of soft materials which enables us to describe 
a unique "rheological fingerprint" of an a priori unknown substance. 

Both the elastic and viscous characteristics of a material can be examined simultaneously 
by imposing an oscillatory shear strain, 7(t) = 7osin(ct;t), which consequently imposes a 
phase-shifted strain-rate 7(^)- Here u is the imposed oscillation frequency, 70 is the maximum 
strain amplitude and t is time. At small strain amplitudes when the response is linear, the 
material is commonly characterized by the viscoelastic moduli G'{uj), G"{uj), as determined 
from the components of the stress in phase with 7(t) and 7(t), respectively. For a purely 
elastic linear solid, the elastic modulus G' is equivalent to the shear modulus G. Similarly, for 
a purely viscous Newtonian fluid with viscosity fi, the loss modulus G" = fiu. However, these 
viscoelastic moduli are not uniquely deflned once the material response becomes nonlinear, 
since higher order harmonics emerge. For convenience the moduli are often determined 
by the coefficients of the first harmonic, G'l and (see Eqn. [T|. These measures of the 
viscoelastic moduli are arbitrary and often fail to capture the rich nonlinearities that appear 
in the raw data signal [S]. 

An example of such rich behavior is shown in the large amplitude oscillatory shear (LAOS) 
results from a wormlike micelle solution in Fig. [T} The periodic stress response a{t;uj,'jQ) 
at steady state is plotted against either •yit) or 7(t), the simultaneous phase-shifted inputs. 
These parametric plots are commonly called Lissajous curves (or more accurately, Bowditch- 
Lissajous curves [21] )■ In this parameter space, a linear viscoelastic response appears as an 
ellipse which is progressively distorted by material nonlinearity. We refer to the a{t) vs. 7(t) 
curves (Fig. [l^) as elastic Lissajous curves to distinguish them from the viscous Lissajous 
curves (Fig. [T|d) which plot a(t) as a function of the shear-rate 7(t). 

The most common method of quantifying LAOS tests is Fourier transform (FT) rheology 
P]. For a sinusoidal strain input 7(t) = 70 sin tut, the stress response can be represented as 
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FIG. 1: Viscoelastic material response from nine experimental oscillatory tests of the micellar so- 
lution. Each trajectory is positioned according to the imposed values {uj, 70}. (a) Elastic Lissajous 
curves: solid lines are total stress cr{t)/amax vs. 7(t)/7o, dashed lines are elastic stress cr'{t)/amax 
vs. 7(t)/7o; (b) Viscous Lissajous curves: solid lines are total stress cr{t)/amax vs. 7(t)/7o, dotted 
lines are viscous stress (t" {t) / a„iax vs. 7(t)/7o • The maximum stress, (Tmax-, is indicated above 
each curve. 



a Fourier series 



a (t; UJ, 7o) = 7o { G"„ {uj, 70) sin nujt 



n odd 



+ G'^ {uj, 'Jo) COS nujt }. 



(1) 



This expression emphasizes elasticity. Viscous scaling results by factoring out 70 = 70a; 
rather than 70, in which case the coefficients are rj'^ = G'^/uj and rj'^ = G'^/u. In the 
linear regime r]' is known as the dynamic viscosity. Only odd-harmonics are included in this 
representation since the stress response is symmetric with respect to shear strain or shear- 
rate, i.e. the material response is the same in both shear directions. Even-harmonic terms 
can be observed in transient responses, secondary flows [ID], or dynamic wall slip [TT], but 
these conditions will not be considered here. Although this FT framework is mathematically 
robust and reduces to the linear viscoelastic framework in the limit of small strains, it lacks 
a physical interpretation of the higher-order coefficients. 

Other methods have been used to quantify viscoelastic nonlinearities [121 1131 [H 115] . 



however these techniques either lack physical interpretation, cannot be calculated uniquely, 
or do not apply simultaneously to both elastic and viscous phenomena. For example, a re- 
cently proposed decomposition [15] uses sets of sine, square, and triangular waves to describe 
a nonlinear response. These basis functions are not orthogonal and thus blur measurement 
of the linear viscoelastic response and qualitative interpretations of the progressive onset of 
nonlinearity. Similarly, the stress decomposition of [H] suffers from non-orthogonality, as 
outlined below. It is therefore desirable to develop a complete and systematic framework 
for quantifying nonlinear viscoelasticity which avoids these ambiguities. 

To illustrate our proposed new framework, we apply it to a wormlike micelle solution and 
a natural biopolymeric hydrogel, gastropod pedal mucus. The wormlike micelle solution 
(prepared as in [16]) consists of cetylpyridinium chloride (CPyCl) and sodium salicylate 
(NaSal) dissolved in brine. The ratio of CPyCl/NaSal is 100 mM/50 mM (3.2 wt%/0.76 
wt%) in a 100 mM (0.56 wt%) NaCl aqueous solution. Pedal mucus samples were collected 
as in [8J. All experiments were performed on a strain-controlled ARES rheometer (TA 
Instruments) equipped with a Peltier plate maintained at T = 22°C, using a solvent trap to 
inhibit evaporation. The micellar solution was tested with a cone-plate geometry (diameter 
D = 50 mm, angle a = 2.3°). The mucus was tested with a plate-plate configuration 
{D = 8 mm and gap h = 550 /im). To eliminate slip, the plate-plate surfaces were covered 
with adhesive-backed waterproof sandpaper, 600 grit (ARC Abrasives Inc.). 

To interpret the data, we extend the method of orthogonal stress decomposition [H], 
which uses symmetry arguments to decompose the generic nonlinear stress response into a 
superposition of an elastic stress (t'{x), where x = 7/70 = sin tut, and viscous stress cr"{y) 
where y = 7/70 = cos tut. Thus, plotting a' vs. x or a" vs. y produces single-valued functions 
in contrast to the closed loops formed by total stress a vs. 7 or a vs. 7 (see Fig. |l|. The 
intra-cjcle elastic and viscous nonlinearities (i.e. nonlinearities within a given steady state 
cycle) are therefore related to the nonlinearity of these functional forms. Cho et al. 
suggest a polynomial regression fit to these lines of elastic and viscous stress. However, the 
material properties represented by these coefficients are not unique since they depend on 
the number of fitting coefficients arbitrarily chosen by the user. 

Instead, we suggest that these curves be represented by the set of Chebyshev polynomials 
of the first kind as they are symmetric, bounded and orthogonal on the finite domain, 
— 1 < X < 1, and can be easily converted to the Fourier coefficients which have dominated 
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the discussion on quantitative LAOS analysis. Using this basis set, the elastic and viscous 
contributions to the measured stress response can be written as 

cr'(x) = 7o ^ e„ {u, 7o) T„(x) (2a) 

n odd 

(^"(y) =io^Vn (t^, 7o) Tn{y) (2b) 

n odd 

where T„(x) is the nth-order Chebyshev polynomial of the first kind. We refer to e„ as the 
elastic Chebyshev coefficients and f„ as the viscous Chebyshev coefficients. 

In the linear regime Eq. [2] recovers the linear viscoelastic result such that ei = G\ and 
^1 = Vi = G'l/uj. We interpret any initial deviation from linearity, i.e. the n = 3 harmonic, 
as follows. A positive contribution of the third-order polynomial T^lx) = Ax^ — 3x results in 
a higher stress at maximum strain, x — f , than represented by the first-order contribution 
alone. Thus, 63 > corresponds to intra-cjcle strain-stiffening of the elastic stress, whereas 
63 < indicates strain-softening. Similarly, a positive value for ^3 represents intra-cycle 
shear-thickening and f 3 < describes shear-thinning. Note that this framework gives a 
direct physical interpretation of the commonly reported third-order Fourier coefficients G3, 
G3 [HI |T7]. As defined in Eq. |2| the Chebyshev coefficients can be related back to the 
Fourier coefficients of Eq. [l] Using the identity T„(cos^) = cos{n9), together with sin 6' = 

n-l 

cos(7r/2 — 9) to show that T„ (sin 6') = sin (nO) (—1) ^ (for n odd), the relationship between 
Chebyshev coefficients in the strain or strain-rate domain and Fourier coefficients in the 
time domain is given by 



n-l 



e„ = GU-1)^ (3a) 

Vn=^=Vn (3b) 

(for n odd). Thus, just as the third-order Chebyshev coefficients have the physical inter- 
pretation described above, so too the third-order Fourier coefficients (with appropriate sign 
corrections for Gg) give physical insight into the nature of deviation from linear viscoelas- 
ticity. When applied to the micelle data of Fig. [T| the Chebyshev coefficients offer physical 
insight to the intra-cycle nonlinearities. For uj = 0.3 rad.s"^ and 70 = {1,3.16,10}, the 
results are 63 = {0.023,0.127,0.235} Pa and V3 = {-0.011,-0.079,-0.288} Pa.s, indicat- 
ing progressive intra-cycle strain-stiffening in the elastic response (63 > 0) combined with 
mira-cycle shear-thinning in the viscous response (^3 < 0). 



This physical interpretation of higher order coefficients motivates several new geomet- 
ric measures for reporting the magnitudes of first-order (linear) viscoelastic moduli in the 
nonlinear regime, to complement the often reported ffist-order Fourier coefficients. These 
additional measures all reduce to G' and G" in the linear regime, but diverge systematically 
when used to analyze a nonlinear signal, offering additional physical insight (beyond that 
captured by the average measures G'l, G'{) into the underlying rheological response. The 
variation in these new measures can be reported as a function of imposed strain amplitude 
7o to indicate the nature of the nonlinearity across different steady state cycles (mter-cycle 
nonlinearities) . 

The nonlinear viscoelastic response of native pedal mucus gel secreted by the terrestrial 
slug Limax maximus provides a striking example of the shortcomings of the conventional FT 
rheology framework [8]. When tested in oscillatory shear, the "elastic modulus" reported 
by the instrument (i.e. the ffist harmonic elastic modulus, G'^) decreases slightly with strain 
amplitude, implying a weak strain-softening. However, examining the raw data in the form 
of Lissajous curves indicates a pronounced local strain-stiffening within a given steady-state 
cycle at sufficiently large stress amplitudes [Hj . A series of strain-controlled oscillatory tests 
at a fixed frequency of = 3 rad.s"^ showing this peculiar nonlinear behavior of pedal 
mucus is given in Fig. |2} The elastic Lissajous curves in Fig. [2]d are elliptical for small 70 
(see inset), but become progressively distorted in the nonlinear regime. At large strains, 
the shear stress is greater than one would expect by projecting the center portion of the 
ellipse, suggesting elastic strain-stiffening which is not captured by the ffist harmonic elastic 
modulus. This behavior is not unique to pedal mucus but appears to be common in soft 
biological materials and can been seen, for example, in data reported by ^18j| for a keratin 
filament network. In both cases intra-cjcle strain-stiffening is readily apparent in the elastic 
Lissajous curves even though the familiar "viscoelastic moduli" do not appear to increase 
as a function of strain amplitude. 

Reporting an "elastic modulus" of such a material as G'l is misleading, as other har- 
monics may also store energy [19]. The ffist-harmonic represents a sine transform, G'l = 
u / (TcyQ) j a(t)'y(t)dt, which is a measure of average elasticity, and is unable to distinctly 
represent the local elastic response of a material at small and large strains. To capture this 
local behavior, we define a set of geometrically-motivated elastic moduli and derive their 
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FIG. 2: Oscillatory strain sweeps of pedal mucus from Limax maximus at a frequency oj = 
3 rad.s^^. (a) Typical rheometer output of the fluid viscoelasticity as parameterized by the first- 
order Fourier moduli, (b) Plotting the raw data from every-other point as a{t) vs. 7(t), with 
graphical representation of elastic moduli G\j,G'^ for 70 = 2.8. 



relation to the Fourier and Chebyshev coefficients. Consider the following: 

da > 



G 1 



M 



= nG'^ = ei - 3e3 + ... 

7=0 n odd 



a 



7 



n-l 
1 2 



(4) 
(5) 



= =61 + 63 + .. 

7=70 n odd 

where G'j^ is the minimum- strain modulus and G'^^ is the large-strain modulus. These mea- 
sures can be visualized graphically as shown in Fig. ^jp by the broken lines. The measures 
are deliberately chosen such that both converge to the linear elastic modulus G' in the limit 
of small strains. We also define a set of dynamic viscosities for reporting the viscous or 
dissipative response. The definitions and relation to the Fourier series of Eq. [T] are 

da 



d'j 



- 

7=0 n odd 



n 

n— 
ijj 



n-l 
-1) — 



V\ - 3^3 + ... 



a 



Vl = 



1 



G" 

= V — = + t^s + 

. . ^-^ OJ 

7=70 n odd 



(6) 
(7) 



The mter-cycle variations of these new measures, as applied to pedal mucus, are reported 
in Fig. [3j Here the minimum-strain modulus softens, G'j^ decreasing with 70, whereas the 
large-strain elasticity first softens then stiffens, G'^ first decreasing then increasing with 70 
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FIG. 3: Oscillatory shear test from Fig.[2| analyzed within the new framework, (a) Elastic moduli: 
minimum-strain and large-strain elastic moduli compared to first harmonic elastic modulus, (b) 
Dynamic viscosities: minimum-rate and large-rate dynamic viscosities compared to first harmonic 
dynamic viscosity, (c) Elastic nonlinearity measures: scaled 3rd order elastic Chebyshev coefficient 
es/ei and strain-stiffening ratio S, both indicate intra-cycle strain-stiffening (d) Viscous nonlin- 
earity measures: scaled 3rd order viscous Chebyshev coefficient v^/vi and shear-thickening ratio 
T, both indicate intra-cycle shear-thinning. 



(see Fig. |3p,). From Fig. [3p it is also apparent that the dissipation at small strain-rates (?7a/) 
increases with increasing 70, whereas the dissipative nature at large strain-rates (represented 
by 7]'^) decreases from one cycle to the next. This rich behavior is obscured by the average 
measures of elastic modulus G'^ and dynamic viscosity •)][ = G"/uj which are commonly 
reported. 

The intra-cycle nonlinearity which distorts the linear viscoelastic ellipse can also be quan- 
tified by comparing these new material measures. Here we define the strain-stiffening ratio 
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as 



G'^ ei + e3 + ... ^ ^ 

where 5 = for a linear elastic response, S > indicates intra-cjcle strain-stiffening, 

and 5* < corresponds to intra-cjcle strain-softening. Users of this framework may also 
choose to compare the moduli using e.g. the ratio G'j^/G'j^ = {1 — S)^^; we chose the former 
for convenience in relating the new measures to higher order Fourier/Chebyshev coefficients, 
and to eliminate potential singularities as G'j^ — > 0. We similarly define the shear-thickening 
ratio as 

y ^ ii - ihi ^ 4^73 + ... 

v'l V1 + V3 + ... 

where T = indicates a single harmonic linear viscous response, T > represents intra- 
cjcle shear-thickening, and T < intra-cjcle shear-thinning. The intra-cjcle nonlinearities 
of pedal mucus are quantified in Fig. ^,d. Both S and 63 are positive at the largest strain 
amplitudes, indicating strain-stiffening, while the nonlinear viscous measures T and are 
both negative at large strain-amplitudes, indicating intra-cjcle shear-thinning. 

In conclusion, the growing interest in biological and other soft materials compels a need for 
consistent, quantitative, low-dimensional descriptions of nonlinear material responses. We 
propose a new comprehensive framework to quantify such nonlinear viscoelastic behavior. 
The scheme provides a physical interpretation of deviations from linear viscoelastic behavior, 
describing elastic and viscous nonlinearities separately, simultaneously, and more thoroughly 
than currently reported measures (summarized by example in Fig. [3]). In addition, the 
method can be easily applied to previously collected data using existing software [22] . These 
new measures may lend insight in the development of constitutive models, and may also serve 
as a more sensitive test for comparing the distinguishing features of different materials. The 
framework is broadly applicable to any complex fluid or soft material which can be tested in 
oscillatory shear, and serves as a complement to the familiar and successful linear viscoelastic 
framework embodied in G'(uj) and G"(uj). 
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